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Abstract: 

Chen and Dohm predicted theoretically in 2004 that the widely believed 
universality principle is violated in the Ising model on the simple cubic lattice 
with more than only six nearest neighbours. Schulte and Drope by Monte 
Carlo simulations found such violation, but not in the predicted direction. 
Selke and Shchur tested the square lattice. Here we check only this univer- 
sality for the susceptibility ratio near the critical point. For this purpose we 
study first the standard Ising model on a simple cubic lattice with six nearest 
neighbours, then with six nearest and twelve next-nearest neighbours, and 
compare the results with the Chen-Dohm lattice of six nearest neighbours 
and only half of the twelve next-nearest neighbours. We do not confirm the 
violation of universality found by Schulte and Drope in the susceptibility 
ratio. 

1. Introduction 

To study the critical phenomena of any system, all systems are divided 
[1] into a small number of universality classes. They are characterized by 
the dimensionality of the space and the number of components of the or- 
der parameter. Within a certain universality class, the universal quantities 
(critical exponents, amplitude ratios, and scaling functions) are independent 
of microscopic details, such as the particular type of interactions or lattice 
structure. 

Once the universal quantities of a universality class are known, the asymp- 
totic critical behavior of very different systems (e.g. fluids and magnets) is 
believed to be known completely provided that only two nonuniversal ampli- 
tudes and the non-universal critical temperature Tc are given. 
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Susceptibility versus |T-Tc|/Tc for L=50(+) and L=1 00(x) for same time =1 0000 
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Figure 1: Susceptibility versus temperature difference for L — 50 and 100 for 
6 nearest neighbours as log-log plot. The upper data correspond to T > Tc, 
the lower to T < Tc. Always 1000 iterations were made. The straight lines 
in Figs. 1,3,4 have the theoretical slope —1.24. 

Here the 3D Ising universality class is considered by studying the suscep- 
tibility of nearest neighbour Ising model (NN model), and the next-nearest 
neighbour Ising model (NNN model) with only six NNN. Here some devia- 
tions from universality were predicted [2] and partially confirmed [3,4]. Also 
with directed interactions problems occur in the Ising model [5,6,7]. Note 
that Chen and Dohm [2] made no prediction on the susceptibility ratio [8] 
and thus are not directly tested in the present paper, which only checks on 
a result of [3] . 

In the following section, we study the susceptibility for the NN model 
and the ratio of the susceptibilities above to below Tc, and then we simulate 
the NNN model for 3D Ising model for different interaction ratios R, all 
with Glauber kinetics. Periodic and helical boundary conditions were used 
throughout. The susceptibihties were calculated through the fluctuations of 
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Ratio of Susceptibility versus |1-T/Tc| for L=50(+) and L=100(x) 
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Figure 2: Ratio of susceptibilities above to below Tc, versus |1 — T/Tc\, for 
the two lattices L = 50 and 100 for 6 nearest neighbours. 

the magnetisation while Ref.3 used the those of the absolute value of the 
magnetisation. 

2. NN Ising Model for 3D 

To test the variation of susceptibility within the 3D Ising universality 
class, we choose first the NN model without external field on the simple 
cubic lattice of L x L x L spins. Each lattice site has 6 nearest neighbours. 
We take two sizes of lattices (L = 50 and 100), where L is the size of lattice, 
and then plot the susceptibihty versus \T — Tc\/Tc double-logarithmically in 
Fig.l. As we see there is little difference in the susceptibility between the 
two lattice sizes. 

Fig. 2 shows the ratio of susceptibilities above and below for both lat- 
tices, L = 50 and 100, at the same distance from T^,, versus this distance. It 
is consistent with the well established value near 5, and thus confirms our 
simulation and analysis methods. 
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L-t : 50-1 00000{+), 1 00-20000(x), 1 00-2000(*), 1 00-1 0000(full sq.), 201 -1 000{empty sq.) 
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L = 40(+), 101 (X), 201 r),301 (empty sq.), 501 (fullsq.), up to 200,000 iterations; Tc{R)=0.6565 
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Figure 3: As Fig.l but with 6 NNN neighbours added to the 6 NN neighbours; 
interaction ratio R — —0.237. The lower part includes larger systems. 



3. NNN Ising Model for 3D 

Also in the case of 6 NN and 12 NNN interactions of equal strength, 
we find a susceptibility ratio near 5 (not shown). Furthermore we study 
the universality of susceptibilities in the zero field NNN Ising model with 
antiferromagnetic anisotropic next-nearest neighbour coupling, added to the 
ferromagnetic NN coupling with a negative ratio R of NNN to NN interaction 
strength. The anisotropic NNN Ising model [2] is established by considering 
only 6 of the 12 next-nearest neighbours being effective for NNN interaction, 
and the other 6 NNN having no interaction. The interacting NNN have the 
position differences ±(1, 1, 0), ±(1, 0, 1), ±(0, 1, 1) on a simple cubic lattice. 
And we do simulations for different sizes of lattices, times and temperatures. 
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L=101 (+,x) and 51 (*,sq.) for t=2000 (+,sq.) and 20000 (x,*) at ratio R = - 0.245, Tc(R) = 0.644 



10 





1 1 1 1 

X 




1 


' 1 . 






X 










>^ 


x~.^ 
X --.^ 










^■^■^ + + 

X, X 










X ": 




... 1 




1 


1 




0.001 




0.01 
|T-Tc|/Tc 


0.1 



Figure 4: As Fig.3 but for R = -0.245, as in [3]. 

and same coupling ratios R — —0.237 and —0.245 as used as the extreme 
cases in [3]. Sizes between 40 and 501 and times between 1000 and 200,000 
are given in the headhncs of Fig.3. The critical temperatures were determined 
by the maxima of the susceptibility, like O.6565 at R + -0.237, in units of Tc 
for the NN model. 

Fig.3 shows the susceptibihty versus the temperature difference to Tc{R) 
as a log-log plot. Since the resulting susceptibility ratio is consistent with the 
one shown for only nearest neighbours and not consistent with the increased 
value found by Schulte and Drope [3], we repeated the simulations at a 
different place using a different computer and larger lattices up to L = 501; 
the results in Figs. 3 and 4 confirm the universality of the susceptibility ratio. 

In general one lattice only was simulated, and the errors thus can best 
be seen by comparing our different symbols at the same temperature, corre- 
sponding to different times and different size; then one also sees the system- 
atic errors. Fig. 3b shows that the data closest to Tc are too much influenced 
by these systematic errors, while the two upper full squares in that Fig. 3b 
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at a relative temperature difference of 0.0015 correspond to L = 501 with 
100,000 (bad) and 200,000 (better) iterations. The critical temperature was 
determined from the maximum of the susceptibility. 

Ref.3 only simulated temperatures 3 % above and below Tc which is in the 
asymptotic regime for the nearest-neighbour case Fig.l but not for the more 
complicated lattice below Tc in Fig. 3. Perhaps that is the reason why they 
got too high susceptibility ratios; our data using a wide range of temperature 
differences show that 1 % difference would have been better and that at this 
smaller difference no reliable deviation from a universal ratio exists. 

4. Conclusion 

The susceptibility ratio, as we see from our simulations for different R = 
0, —0.237, —0.245, keeps the universality of the 3D Ising model above and 
below the critical temperature. For 6 nearest neighbours, for 6+12 nearest 
next-nearest neighbours and for 6+6 nearest neighbours, there is no reliable 
difference for the susceptibility ratio, contradicting [3] but compatible with 
the standard ratio of 4.75 [9]- We did not attempt to simulate the more 
subtle universality questions connected with Binder cumulants [2,4]. 

We thank W. Selke for comments on the manuscript. DS thanks Univer- 
sidade Federal Fluminense in Brazil for its hospitahty during the time of his 
simulations. 
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